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Abstract 

We study the behavior of cold atoms trapped in optical double well potentials in 
presence of noise, either generated by an external environment or by the trap itself. 
We show that quite in general the noise can induce a current between the two wells even 
when tunneling between the two wells is highly suppressed (Mott insulator phase). An 
engineered environment could provide a test ground for the behavior of such a current. 



1 Introduction 

Cold atoms (typically Bose-Einstein condensates) in optical lattices provide a very man- 
ageable and versatile tools for experimental tests in many-body quantum physics and as 
such have recently attracted a lot of interest In a suitable approximation, the dynamics of 
bosonic atoms in this kind of traps can be described by a Bose-Hubbard model [2], which, 
depending on the value of the ratio between the hopping amplitude and the on site repul- 
sive boson-boson interaction, predicts a superfluid and a Mott insulator phase. Much of the 
experimental work on such systems have been dedicated to the study of this phase transition. 

These systems are usually treated as isolated. However this is only an approximation, and 
there are instances when the presence of a weakly coupled external environment can not be 
neglected. Typically, this happens when the optical lattice is immersed in a thermal bath, 
either external or formed by the fraction of cold atoms that are not in the condensed phase, 
or when the lattice itself is subjected to stochastic noise. In these cases, the condensate gas 
in the lattice need to be treated as an open quantum system, and its dynamics is altered by 
effects of noise and dissipation [21 IH El El El E] • 

In the following, we shall study such effects in the particular case of a lattice with just 
two sites, where the cloud of condensed bosonic atoms is confined in a double-well potential 



^The literature on this topic is enormous; for a recent review see [T] and the long list of references therein. 
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(extension to the multiple sites case will be reported elsewhere). In the ideal, isolated 
situation, the Mott insulator phase corresponds to a blocked regime where the occupation 
numbers of the two wells is fixed and no atom is allowed to cross the barrier passing from a 
well to the other. In particular, no current across the trap barrier is observed. 

On the other hand, when the trap can not be considered isolated, the two wells can be 
effectively connected through the indirect action of the environment, whose presence may 
then break the insulating regime. This is in fact what it is found by analyzing in detail the 
time evolution of the current operator: as soon as the interaction with the environment is 
switched on, an environment induced current is generated via a purely noisy mechanism. 
Through an engineered environment, e.g. by injecting and modulating a stochastic noise in 
the double well potential, this current can be in principle experimentally analyzed. 

2 Double- well trap 

As explained above, we shall study the behavior of a system of condensed bosonic atoms 
trapped in a double well potential, weakly coupled to an external environment. In a suitable 
approximation, the dynamics of the system alone {i.e. in absence of the environment) is 
very well described by a Bose-Hubbard type hamiltonian |2j, that in the present case takes 
the simple form 

Hs = -T{a[a2 + alai) + U{nl + n\) + eiUi + e2n2 , (1) 

where the operators Oj, a\ , i = 1,2 are annihilation and creation operators for atoms in 
the two wells, obeying the standard oscillator algebra [oj, a]] = 6ij, with = a|aj, the 
number operators. The various contributions in the hamiltonian ([T]) are characterized by 
their physical meaning: the first is an hopping term dependent on the tunneling amplitude T, 
while the remaining are on-site energy terms; the one proportional to the coupling constant 
U is quadratic in the number operator and is due to the boson-boson repulsive contact 
interaction, while the second one is due to the trapping potential and may be different for 
each well. The parameter ei represents the energy of the bottom of the i-th well and therefore 
describes its depth. 

As well known [1], the hamiltonian ([1]) describes the quantum phase transition between a 
superfluid and insulator phases, whose order parameter is given by the ratio T/U. We are 
interested in studying the regime where the atomic cloud is well separated in the two wells; 
this is the case when the hopping energy is very small compared to the on-site one. Neglecting 
the shifts Ei which are of the same order for both the wells, when an atom hops, it loses 
an energy T and gains instead an energy U; when T/U <^ 1, this situation is energetically 
suppressed and, at equilibrium, the ground state of the system is given by a Fock state with 
equal number of bosons per well, |A^, A^) = |A^) ® |A^). This is the typical situation of a 
Mott insulator, where no current is flowing between the two wells. 

In order to give a precise meaning to this statement, one needs to introduce a suitable 
current operator. Following P [ID] , we shall first consider the barycenter of the system of 
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atoms in the two wells; in suitable units, it is simply given by 

{a\ai - a\a2) 



2N 



(2) 



The velocity with which the barycenter moves can be obtained taking the time derivative of 
this operator: 

— = i[Hs, Z] = ^(«1«2 - 4«i) • (3) 

The current operator describes the flow of atoms between the two wells and thus it should 
be proportional to this velocity operator; in the chosen units, it can be simply taken to be 

J = i{a\a2 — alai) . (4) 

In the ground state \N, N), this operator have vanishing mean value, meaning that no 
current flows across the barrier when the wells are equally filled. Moreover it is easy to check 
that also its time derivative vanishes: 

^^{J)^ = ^^[HBH, J])^ = ^- (5) 
As we shall see, this conclusion is in general no longer true in presence of an environment. 



3 Noise induced current 

When the interaction with the environment E is switched on, the hamiltonian describing the 
evolution of the total system, condensed gas in the trap plus environment, can be decomposed 
as 

HsE = Hs + He + XHj , (6) 

where Hs, driving the free motion of the atoms in the double well trap, is as in ([T]), He de- 
scribes the evolution of the environment alone, while Hj takes care of the interaction between 
atoms and environment, A being a small coupling constant, A <^ 1. For a weakly coupled 
environment, the form of Hj can be taken to be bilinear in the system and environment 
variables, 

4 

Hi = J2yi®Bi, (7) 
1=1 

where Bi are hermitian operators describing environment observables, while the system op- 
erators Vi are given by the hermitian combinations of the oscillator variables: 

Vi = {(fli + a{), i{ai - a\), (02 + 4)^ '^('^2 - 4)} • (§) 

The state of the composite system S + E, described in general by a density matrix psE, 
evolves unitarily with the hamiltonian However, the dynamics of the density matrix 
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p = Tte[pse] describing the state of the condensed atoms in the trap and obtained by 
averaging over the environment degrees of freedom, is in general comphcated, containing non- 
hnearities, secular terms and memory effects [51 HJ El El El E] • When system and environment 
are initially uncorrelated, no secular terms or non-linearities arise; this is a rather common 
situations in actual experimental realizations, so that we shall henceforth assume that at the 
initial time t = 0, the state of the composite system be in factorized form, Pse{0) = p{0)'^Pe, 
Pe being the reference equilibrium state of the environment, e.g. a thermal state pe — e~^^'^ . 

Memory effects are typically short time phenomena, due to the initial interaction between 
subsystem and environment. We are not interested in the details of such phenomena; rather, 
we would like to study the effective dynamics of the subsystem alone on long time scales, for 
which memory effects have already died out. Operatively, since the noise and dissipation due 
to the environment start to become relevant on time scales of order A~^, this is obtained by 
passing to a suitable rescaled time variable, t — »• t/A^ and then letting A — 0. This procedure 
is known as "Markovian approximation" and allows disregarding the initial memory effects. 
In order to be physically applicable, this approximation requires an a priori unambiguous 
separation between subsystem S and environment E. In the weak coupling regime, this is 
typically achieved when the ratio te/tr, between the typical decay time of correlations in 
the environment and the characteristic relaxation time of the subsystem immersed in the 
environment, is small. 

Two different Markovian limits have been discussed in the literature, depending on the 
way of rendering the ratio teItr small. In the weak coupling limit [HI [T2], the correlations 
in the environment 

G,j{t) = {B,{t)B,) = Tr [B,{t) B, pe] , Bi{t) = e'*^^5, e"'*^^ , (9) 
are non-singular and assumed to decay fast enough, 

POO 

/ dt\G ij {t)\{l + ty < oo , (10) 
Jo 

with e a positive constant. For typical environments, like heat baths, one finds an exponential 
behavior: 

Gi,it)-G,j e-^^\ Pe = 1/te, (11) 

with Gij a constant hermitian matrix. In the singular coupling limit [13] instead, it is the 
decay time of correlations in the environment that becomes small, te — »■ 0; concretely, this 
is obtained by letting the correlation functions tend to a Dirac-delta: 

G,,it)^G,,6{t) . (12) 

This situation is encountered in a bath with infinitely large temperature or in the case in 
which the subsystem is subjected to stochastic noise; in the latter case the positive constants 
Gij form a symmetric matrix. 
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In both cases, the master equation describing the evolution of the subsystem density 
matrix p{t) takes the form 

^ = Up{t)] ^ -i[Hs + H^'\ pit)] + D[p(t)] . (13) 

The environment modifies the standard Liouville-von Neuman equation for p by an effective 
correction if*^^) to the free hamiltonian Hs and a pure dissipative contribution D which 
cannot be written in hamiltonian form. Explicitly, these new contributions take the standard 
Kossakowski-Lindblad expression [TU [151 [E] : 

4 



Ld ij=l 



V^^ico)pV,ico)-Uv,{cu)V^H^),P 



(14) 



uj ij=l 

where the formal sum over ut represents the sum over all possible energy differences within the 
spectrum of the subsystem hamiltonian Hs- The two hermitian matrices, Cij{uj) and Sij{uj), 
are the Fourier and Hilbert transform of the environment correlation functions, respectively, 
{V denotes principal value), 

/+00 
dte-*"*G,,(t) , (16) 
-co 

s.,ico) = ^V r^dw^. (17) 

They embody the effects of noise and dissipation induced by the environment. Notice that 
the Kossakowski matrix Cij{u!) turns out to be positive, a key requirement in order to guar- 
antee the positivity of p{t) for any time, and thus the physical consistency of the evolution 
generated by f|T3|) [3l, (8]. Finally, the Kraus operators Vi{u) are suitable sandwiches of the 
system operators \^ of ([8]) between eigenprojectors of the system hamiltonian Hs (see the 
following for explicit expressions). 

By exponentiation, the operator L[-] generates a family of completely positive and trace- 
preserving maps t — > 7t = e*'^ which constitutes a so-called quantum dynamical semigroup 
[21 [1], for which: '^t°ls = It+s-, ^, -s > 0, thus encoding the intrinsic irreversibility of the 
dynamics. 

As discussed before, we are interested in analyzing the behavior of the current operator 
^ in presence of the environment. Its evolution in time is generated by a master equation 
which is obtained from f[T3l) by duality: 

^ = L*[J(t)] = i[Hs + H^'\ J{t)] + W[J{t)] , (18) 

where the operator L*[J] is defined via the duality relation Tr(L[p] J) = Tr(pL*[J]), and 
similarly for D*. We shall now study in detail the behaviour of the current operator J{t) in 
presence of noise as described by the two regimes of weak and singular coupling limit. 
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3.1 Weak coupling limit 



As discussed before, we are considering the situation in which, prior to the couphng with 
the environment, the double well trap is in the Mott insulator phase, in which tunnelling 
as induced by the first term in the system hamiltonian ([T]) is suppressed, i.e. T ^ U. The 
noise contributions generated by the presence of the environment are second order in the 
interaction hamiltonian Hj and therefore turn out to be of order (c/. the expressions in 
(fT^ . (fTSjl ): in evaluating the coefficients (fT6|) . (fT7|l and the energy differences uj in the sums 
(fT^ . (fTSl) . one can then safely neglect tunneling contributions, being subdominant, of order 
T. In this approximation, the spectrum of the system hamiltonian ([1]) results quadratic 
in the the occupation numbers of the two wells, while the corresponding energy differences 
can be labelled by an integer n, 



OJr, 



e + U + 2Un . 



£2 



(19) 



e; we 



For simplicity, we are focusing on a trap formed by wells with equal depth, ei 
will comment on the case ei 7^ £2 later on. 

Explicit evaluation allows casting the dissipative contributions to the dual master equation 
(fT8|) in the form: 



r[.j] 



n ij 



v,(cu„) jv;(cu. 



i{Vi(cu„)V](^„),j} 



the Kraus operators are now given by 



1, l®P„a2, 1 



(20) 
(21) 

(22) 



where P„ are projector operators on single well Fock states with occupation number n, while 
the Kossakowski matrix takes the form 

//iii(u;„) 




h{uJn) 





h22i-UJn 



h 
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h24,{—UJn 





(23) 





h2i{~^n 

The entries of this matrix are linear combinations of the Fourier transforms of the correlation 
functions given in (fT6|) : omitting for simplicity the Un dependence, one has: 



(24) 



hii 


= Cii + C22 + i(c21 


- C12) 


hi3 


= Ci3 + C24 + i{c23 


- Cm) 


h22 


= Cll + C22 + i{ci2 


- C21) 


h24 


= Ci3 + C24 + i{cu 


- C23) 




= C33 + C44 + i(c43 


- C34) 


/l44 


= C33 + C44 + i{c34 


- C43) 
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The matrix kij has a similar structure in terms of the Hilbert transform of the environment 
correlation functions; it will not be relevant for the considerations that follow and thus, for 
simplicity, its explicit expression is omitted. 

Assume now the system be prepared in the Mott insulator phase, with initial state pN = 
\N, N){N, N\; as already observed, the average of the current operator (jll) in this state 
vanish: {J{0))n = 0. In order to see whether the presence of the environment alters this 
situation, it is sufficient to study the time derivative of the average {Jit))^ at t = as given 
by the evolution equation (ITSl) with (l20l) and (pTI) . The explicit computation gives: 

d{J{t))N 
dt 

since all hamiltonian contributions vanish, with 



= (D*[J(0)])^ , (25) 





2\' 



{N+lf'^m{hsi{ujN))+N^':sm{h24i-uJN^i)) • (26) 



Therefore, as soon as the interaction with the environment is switched on, a current starts 
flowing between the two wells; it has a purely noisy origin and its magnitude is determined 
by the elements of the Kossakowski matrix (1231) . i.e. by the correlation functions in the 
environment. In the case of an heat bath with correlation functions of exponential type as 
in (llip . one flnds that, for small times and large enough N, the current behaves as: 

= '4p^s^e{G,,-G,,)t (27) 
dt n^ + u^j ^ 

^ ^^3?e(Gi4-G23) t. (28) 

The intensity of the current clearly depends on the strength Gij of the environment corre- 
lation functions. In particular, the current vanishes if the coupling of the two wells to the 
environment is realized through the same operators, i.e. if the environment acts exactly in 
the same way on the two wells; indeed, in this case the r.h.s. of fl2S]) is zero since Bi = 
and B2 = -B4, and therefore, being Gij hermitian, 3?e(Gi4) = 3?e(G'23)- 

Finally, let us briefly consider the situation in which the two wells have initially different 
depth, i.e. ei 7^ £2 in the starting hamiltonian ([T]). The weak coupling limit procedure 
involves an average over the microscopic fast motion of the system as generated by ([1]); this 
is perfectly justifled, since we are focusing on the dynamics over rescaled, long time intervals 
(this is usually referred to as "rotating wave approximation" O Hj). When ei 7^ 62, this 
average procedure gives rise to a simplifled Kossakowski matrix, where only the diagonal 
elements are non vanishing. In other terms, everything goes as if the two wells interact 
independently with the environment and no correlations between them could be created 
through noisy effects, or, more precisely, these correlations average to zero if they are observed 
on time intervals rescaled by A^. As a consequence, since the environment can not connect 
the two wells, the system stays in the insulator phase also after coupling to the environment 
and the average value of the current operator remains zero. 
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3.2 Singular coupling limit 



As mentioned before, this Markovian approximation is appropriate when our double well 
trap is immersed in an environment whose correlation functions decay very fast, as happens 
for instance in a heat bath at infinitely large temperature, or when the trap is subjected to 
stochastic noise. In practice, one assumes the fluctuations in the bath to be delta-correlated 
as given in f[T^ . so that the dissipative contributions to the master equation f[T^ become 
independent from the system hamiltonian ([1]) [13]. The sum over the spectral parameter uj 
disappears from the expressions of D and if*^^-* in (fT4|) and f|T5l) . while the Kossakowski matrix 
becomes constant, Cij = Gij, and Kraus operators Vi as in ([8]). Notice that in general all the 
entries of Cij are now non vanishing, so that the structure of the dissipative contributions in 
( IT3ll is in general richer than in the weak coupling limit 

The dual equation, generating the dissipative dynamics for the current operator, takes 
again the form f[T5]) . where the dissipative contributions can be cast in the form similar to 
(I2D]) and (JSP); explicitly: 



v.Jvj--{v.vj, J 



(29) 



H'-'^ = X'J2h,ViV] , (30) 



where now one simply have 

V = ^fli, a\, 02, . 

The entries of the new Kossakowski matrix are linear combinations of the old one 
Cij, with relation analogous to those in fl2Sl) : similar results hold also for the hamiltonian 
contributions Sij. 

In order to see whether also in this case a current can be generated by purely dissipative 
effects, one first prepares the double well trap in the Mott insulator phase, described by 
the equilibrium state = \N, N){N, N\; then, one switches on the interaction with the 
environment and looks again at the behavior of the average {J{t))N for small times. Here 
again only the dissipative part D* contributes so that the result fl25l) is still valid, while the 
explicit evaluation gives: 

{J{t))M ^ ^^d^^h. t = 2X^ [(2iV + 1) Qm{G3i + G42) + - G23)] t . (31) 

Also in the case of delta-correlated environments, a current is initially generated by a purely 
noisy mechanism; its intensity depends again on the entries of the Kossakowski matrix, i.e. 
on the strength of the correlations in the bath. 



^ This is a consequence of the fact that no "rotating wave approximation" is necessary in the singular 
couphng hmit. 
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From the expression in (13T!) . one sees that as in the case of the weak couphng hmit, the 
current vanishes when the environment couples in the same way to the two wells (as before 
3?e((j'i4) = 3?e(G23), while Gm = Gn, 6*42 = G22 are real, again by the hermiticity of Gij). 
On the other hand, the mechanism generating the noisy current contribution is independent 
from the details of the microscopic system dynamics given by the hamiltonian ([1]); this 
implies that a delta-correlated noise can in principle generate a current even for 61 7^ 62, i.e. 
when the two wells are initially at different depths. 

4 Outlook 

An optical double-well trap filled with condensed bosonic atoms immersed in an external, 
weakly coupled environment represents an example of an open quantum system. The pres- 
ence of the environment generates noise and dissipation, and as a consequence, the evolution 
in time of the atoms in the trap is no longer reversible. Rather, it is described by a quantum 
dynamical semigroup, generated by a master equation of the Kossakowski-Lindblad form 

(USD, with ([HD, (US]). 

As discussed above, the detailed structure of such evolution depends crucially on the 
characteristic properties of the environment, specifically on the behaviour of its two-point 
correlation functions. Nevertheless, as a general result, the dissipative dynamics allows for 
an indirect coupling between the atoms in separate wells, through the noisy action of the 
environment. Therefore, even if the trap is initially prepared in the Mott insulator phase, 
where tunneling is not allowed, the subsequent coupling to the environment may allow a 
transitions between the two wells. 

This possibility has been studied by analyzing the behaviour of a suitably defined current 
operator, signaling an unbalanced flow between the two wells. As soon as the interaction 
with the environment is switched on, the current starts developing a non-vanishing, purely 
dissipative contribution to its average. This contribution can be isolated by preparing the 
trap in the insulator phase, where the current has a vanishing expectation value. The 
interaction with the environment will then have the net result of effectively "breaking" the 
insulator regime. 

This effect may be experimentally ascertained through the use of an engineered environ- 
ment. The idea is to inject stochastic noise into the optical trap, operation that, in the 
formalism developed above, can be described by adding to the system hamiltonian ([1]) an 
external Gaussian stochastic potential. One can show [131 E] that the average over the 
noise produces effects similar to those obtained through the couphng with a delta-correlated 
environment, and thus, as in fl?Il) a non vanishing dissipative contribution to the current 
operator. By releasing the trap potential right after the source of stochastic noise has been 
switched on, one should then be able to measure an unbalance in the barycenter position 
of the ballistic expanding cloud of condensed atoms, a clear signal of flow of atoms between 
the two wells. 
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